Long-lived mesoscopic entanglement outside the Lamb-Dicke regime 
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We create entangled states of the spin and motion of a single 40 Ca + ion in a linear ion trap. The 
motional part consists of coherent states of large separation and long coherence time. The states are 
created by driving the motion using counterpropagating laser beams. We theoretically study and 
experimentally observe the behaviour outside the Lamb-Dicke regime, where the trajectory in phase 
space is modified and the coherent states become squeezed. We directly observe the modification of 
the return time of the trajectory, and infer the squeezing. The mesoscopic entanglement is observed 
up to Aa = 5.1 with coherence time 170 /xs and mean phonon excitation n — 16. 



A two-state system interacting with a quantum har- 
monic oscillator has for a long time played a fundamen- 
tal role in quantum optics 1], and more recently has 
attracted interest in the context of mesoscopic quantum 
physics, micromechanical oscillators [2, 3], and quantum 
information 0, Hj]. The states of motion of a quantum 
oscillator which most closely resemble classical states of 
motion are the Glauber coherent states [l|. It has been 
pointed out that a superposition of such states, with a 
large difference between their coherent state parameters 
a, is a state of affairs comparable to the Schrodinger 
Cat thought-experiment [H, @|. This permits an inves- 
tigation of mesoscopic quantum physics using this sys- 
tem. The conditions of the thought-experiment are most 
closely matched when the superposition involves an en- 
tanglement between the large system (the oscillator) and 
a smaller (e.g. two-state) system, viz: 



|V,) = — (|a T ) |T) +e^K)|j» 



(1) 



where Aa 2 = — aj 2 is large, |f) , |J.) are the states of 
the two-state system, and the interference phase ip must 
be stable and under control in the experiment (as must 
the values of ot-\,i, including their relative phase). 

The coherent states are mesoscopic in that their en- 
ergy is S = | a | 2 3> 1 units of the fundamental exci- 
tation energy (the energy gap between the ground and 
first excited states) and the separation x s of the mo- 
tional wavepackets is greater than their individual size 
xq by the ratio x s /xq = 2Aa 1. The size of the 
Hilbert space required to express the motional state is 
approximately log 2 n qubits, and in the case of a state 
such as (TTJ there is entanglement with another degree 
of freedom. The phase coherence time, T2, is sensitive 
to the separation [7j. These measures are summarized by 
the list {n, Aa,x s ,T2}. 

States of the type (p} have been realized in the inter- 



nal state and motion of single trapped ions [8|, [9|, 1C|, 
and in an atom interacting with a cavity- field [111 1 12 1 . 
For experiments where the coherence of the two parts of 
the state was observed [l3j], the reported parameter values 
were {n,Aa,x s ,T 2 } = {8.8, 2.97, 42nm, O(10^s)} d; 



{12, 5.2, 73nm, 6^s|Ji; {1, 2, 14 nm, ~ 0.5ms} 0; 
{9.5, 1.8, — , 90 /is} HJ 

We present experiments in which the cat state is real- 
ized with large values of both the size and the coherence 
time together, and we describe and demonstrate a quali- 
tatively new behaviour which appears outside the Lamb- 
Dicke regime (LDR). We have generated cat states of 
the spin and motion of a single trapped 40 Ca + ion with 
{n, Aa, x s ,T 2 } — {16, 5.1, 170 nm, ~ 170 /is}, observing 
their coherence using an interference effect. The Hilbert 
space dimension is approximately 5 qubits. In our exper- 
iments the driving of the motion goes outside the LDR: 
that is, the force on the ion is not independent of its 
position, in contrast to previous work. This results in a 
dramatic modification of the trajectory in phase space, 



and also squeezing of the coherent state [14| . This means 
that we achieve cat states of the traditional type, and 
also infer production of states in which the entangled and 
subsequently interfering components are not standard co- 
herent states, but squeezed states, with a squeezing pa- 
rameter (ratio of principal axes of the Wigner function) 
~ 3. 

Our system consists of a single spin-half particle in a 
harmonic potential, subject to a "walking wave" of light 
formed by counterpropagating laser beams in a stand- 
ing wave configuration with a frequency difference ap- 
plied between the two beams. The walking wave pro- 
vides a spin-dependent force on the particle. The inter- 
action Hamiltonian is Hj = + ^2 m V m \m) (m\ where 
m =1,1, 



V m = Ml cos(kx — Lut + <f) m ) 



(2) 



and H w = fiA ff (|T) (T| - II) (ID/2. V m is a light shift 
from far-off-resonant single-photon excitation; H T is a 
light shift from off-resonant Raman excitation of spin- 
flip transitions. The latter has no effect on the motion, 
but causes the spin state to precess. k and uj are the 
wavevector along the x axis and the frequency of the 
walking wave respectively. 

The position-dependence of V m results in a spin- 
dependent force f m {x,t) = —dV m /dx. The clas- 
sical equation of motion is 2Mu>oXo(da/dt) = 



2 



iexp{ia>ot)f m (x,t) where a — exp(iujQt)(x + 
ip/Mu )/2xo, and xq is a length scale. M is the 
mass of the ion and loq its natural oscillation frequency 
in the trap. If we take Xq = (h/2AIuj ) 1 / 2 then a 
corresponds exactly to the coherent state parameter in 
the quantum treatment. 

We consider motional states \a) starting at or near 
a = 0. For small | (kx) | we have the LDR, where 
the force f m (x,t) ~ Mlksin(u>t — (f> m ) is independent 
of position. In this case an analytical solution of the 
time-dependent Schrodinger equation (TDSE) is possible 
0, H, Ell ■ The quantum state is merely displaced along 
a trajectory a(i) exactly matching the classical predic- 
tion. For | S | <C where 5 = u — ojq, a(t) describes 
a circle of diameter rfil/8, where r\ = kxQ is the Lamb- 
Dicke parameter. The motion around the circle is at a 
constant rate, completing a loop at integer multiples of 
2ir/5. The quantum state picks up a phase proportional 
to the area of the loop, plus a contribution ±A w i/2, and 
an oscillating phase from J V m dt. This oscillating phase 
scales as 1 / 77; it is important when (kx) <C 1 but is small 
outside the LDR when (kx) ~ 1 and we will ignore it 
hereafter. 

Outside the LDR, the classical dynamics cannot be 
solved analytically. We have studied this by numerical 
analysis: a numerical integration of the classical equation 
of motion, and an approximate numerical solution of the 
TDSE. The latter included up to n max = 100 harmonic 
oscillator levels and terms in all orders of rj for the carrier 
and first three motional sidebands in Hj. 

The position dependence of the force causes the tra- 
jectory to be non-circular and the motional state is 
squeezed, that is, the wavepacket is narrowed in one di- 
rection in phase space and broadened in another. For 
modest values of (kx) we find that the squeezing is 
negligible and the quantum wavepacket simply follows 
the modified classical trajectory. For larger values the 
wavepacket is squeezed, and for larger values still, or 
longer times (e.g. after more than one loop in phase 
space) the wavepacket changes shape dramatically. Fig- 
ure 1 shows an illustrative example. The departure from 
a circular trajectory and the squeezing are clearly seen. 
Note that a(r), where r is the length of time that the 
force is applied, still returns to the origin, but at a time t r 
earlier than the value 2ir/5 for a state which stays within 
the LDR. 

We find a(r) begins to differ clearly from a circle when 
(kx max ) = 2?7|Q! max | > 1. Each loop is shaped like a 
tear-drop and t r is reduced. Fortunately, there are two 
features of the motion which simplify the interpretation 
of our experiments. First, in the early stages, it is still 
within the LDR so the initial behaviour both in terms of 
amplitude and phase is accurately described by the sim- 
ple analytical treatment outlined above. Second, when 
the amplitude of the motion drops for the first time back 
into the LDR the analytical treatment again gives a good 
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FIG. 1: (color online) Trajectories in phase space a(r) (dot- 
ted and dashed lines) and motional Wigner function (contour 
lines at 0,1,2 standard deviations) for the parameter values 
lu /2iy = 536 kHz, rj = 0.244, Q/27T = 93 kHz, 5/2tv = 3.4 kHz 
and varying forcing time r. The two trajectories are shown 
for the different spin states for & m = 1.41, up to the time 
2tt/S — 286 /is. & w = |0| — (pi I is the phase angle between the 
forces on the two spin states. The return time is t r = 192 fis. 
For clarity, the Wigner function is shown for just one of the 
trajectories at three example times, r = 0, 96, 286 /is. The 
squeezing is ~ 3 at t r /2. In the experiments, superpositions 
of motion along both members of such pairs of trajectories 
are created. 



representation of the behaviour, if one takes into account 
the difference between the actual and LDR return times. 

We experimentally investigated the behaviour using a 
single 40 Ca + ion in a linear ion trap 17 1. The two-state 
system is the spin- 1/2 ground state of the ion, and the 
potential @ is realized by the light shift when the ion is 
illuminated by a laser walking wave far (30 GHz) detuned 
from the A = 397 nm S 1 / 2 ~Pi/2 transition. A quantiza- 
tion axis is defined by a 1.4 Gauss magnetic field B at 57° 
to the trap axis x. These axes are horizontal. A 60° pair 
of laser beams forms the walking wave along x, with dif- 
ference wavevector k = 2tt/X. Their difference frequency 
uj is generated with 1 Hz precision by acousto-optic mod- 
ulation. One beam is horizontally polarized, the other 
is close to linear at 69° to the vertical. The resulting 
light field has three components: a <r + polarized walking 
wave, a a~ polarized walking wave, and a predominantly 
7r polarized travelling wave. The transition in the ion is 
J = \ — > |, so the rr + (c - ) light interacts with ||) (|t)), 
giving rise to V±, (Vj-) respectively. 

The ion is first prepared in ||) |no) where |no) is a ther- 
mal motional state close to the ground state with mean 
motional state occupation number no fl8| . A sequence of 
laser pulses is then applied, and finally the spin state is 
measured by selective shelving followed by fluorescence 
detection, see 0. This is repeated 500 times to accu- 
mulate statistics, then a parameter value is changed and 
the sequence repeated. 

Initial experiments were carried out by Ramsey inter- 
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FIG. 2: Observed fringe amplitudes as a function of r, nor- 
malized to the value at r = 0. Each point is obtained from 
a sin fit to a scan over <j>. The lines are fitted curves using 
0, with parameter values given in table 1. Dashed horizon- 
tal lines separate data sets taken at different trap frequencies. 
Dashed vertical lines are drawn at r = 2n/6. 



ferometry, with the oscillating force pulse W applied in 
the gap. The n/2 pulses were Raman transitions driven 
by the walking wave, with tunable relative azimuthal 
phase <j). The data reported here were obtained using 
a spin-echo sequence, to eliminate slow phase noise. The 
first it/2 pulse evolves the spin to (|T) + \{))/\/2. After 
the W pulse, of duration r, the system is (up to a global 
phase) in the state ([1]), with ip = A„t. In the LDR, 
\oL\ t i) are coherent states, and outside this limit they are 
squeezed or more general states, centred in phase space 
close to a-|-,x- The observed signal after the final pulse is 



P(t)= 1-Re 



i(0-A„.r) 



/2- 



(3) 



We determine (a-fla;) by observing -P(T) as a function 
of <j) and r. For each value of r we accumulate the in- 
terference fringe pattern as a function of </>, and fit it 
with a sinusoid. To factor out the effect of magnetic 
field noise, we normalize the observed fringe amplitude 
by comparison to that obtained in a control experiment, 
having exactly the same timing but with no W pulse. 
The amplitude of the control experiment fringes drop to 
typically 40% for a 300/is spin-echo time. 

The observed amplitude A and phase offset </>o of the 
fringes are shown as functions of r for various data sets 
in figures 2,3. This information enables us to infer the 
evolution of the system. The analysis is simplified by 
the fact that critical data exist only where the motion is 
within the LDR; outside this region the fringe amplitude 
is close to zero. To fit A we can therefore use a LDR ex- 
pression pi 0] , modified to take account of the reduction 
in return time: 
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FIG. 3: Observed phase offset of the fringes as a function of 
t, successive data sets are offset by 20 radians. The lines are 
fitted curves using (|6]), with parameter values given in table 
f. 



Here 7 is mainly a decay caused by decoherence effects 0], 
but also includes a contribution due to squeezing. The 
reduction in fringe visibility due to squeezing ranges from 
less than 1% for data set 1 to approximately 9% for data 
set 3. The return time t r constrains the global trajectory, 
and allows our most direct observations of non-Lamb- 
Dicke behaviour, in that we find in general t r is signifi- 
cantly less than 2ir/5. The parameter D is related to ao, 
the maximum a which would occur for the same force if 
the motion were entirely within the LDR. We have 



D = Ra \ / 2n + 1 sin($ u ,/2), 



(5) 



where the two trajectories are separated by the angle 
&w = I^T — 4*1 1 an d -ft = t r 5/2ir is the fractional reduction 
in return time. 

The phase is fitted by a similarly modified LDR ex- 
pression 



(t) = (const) + A^t + B 2 sin 2 (7rr/^), 



(G) 



A(t) = e~ 1T exp(-2D 2 sin 2 (irT/t r )). 



(4) 



where B 2 = R 2 al sin($ w ). 

With the polarization angles in the experiments, $ w = 
1.41(5)rad. The amplitude and phase analyses thus each 
give values of Ra and t r . We are then able, using the 
results of our simulations, to determine R (and hence 
values for the detuning and ao), a m ax and Aa max . As a 
check on the validity of our interpretation we can com- 
pare the results obtained for ao and d/2ir with those 
expected on the basis of our knowledge of the laser field. 
The detuning <5/27r is known to ±0.5 kHz. Two pieces of 
information quantify the light intensity: the Rabi flop- 
ping rate Q c when spin-flip ('carrier') transitions are res- 
onantly driven, and the light shift A T (deduced from ([6|)) 
which comes from off-resonant excitation of these same 
transitions during the W pulse. The laser polarization 
is set to be (very close to) linear at the ion for either 
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5.1 


4 


1.50 


91 


3.5 


2.03 
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2.72 
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0.02 0.245 
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TABLE I: Experimental parameters and results. Column 1 gives the data set number. The next 5 columns give values extracted 
directly from the fringe data by curve fitting. The rest of the table gives further raw information and derived quantities, n is 
known from the trap frequency. fi c is obtained by Rabi flopping on the carrier (value a) and from the fitted A-n- (value b). The 
detuning 8 is set experimentally to within 0.5kHz for a given data set (c) and can be evaluated also (via the TDSE) from the 
data analysis (d). The same process leads to a value of ao (f) which can be compared with that deduced from the parameters 
of the light field (e). Finally, we infer a ma x, the maximum excursion in phase space, and Aa max , the maximum distance in 
phase space between the two spin components, from the TDSE. We estimate b,d,f, a max and Aa max have 5% accuracy; the 
consistency check a,c,e combines Rabi flopping, relative power and polarization measurements and is accurate to ~ 15%. 



beam acting alone, by nulling any differential light shift 
observed in Ramsey experiments. 

For all the data there is reasonable agreement between 
observations and predictions. Sets 1-3 were taken on the 
same day under particularly stable conditions, and have 
very good overall consistency. The largest a max and Aa 
(4.0 and 5.1) were obtained with set 3. In particular, we 
note the large reduction in return time (R = 0.67). The 
results of the numerical solution of the TDSE for this 
case are shown in figure 1. 

The observed coherence is not perfect, owing mainly 
to magnetic field noise, photon scattering, and heating. 
The effect of the first of these is small when we use the 
spin-echo sequence. Let a < 1 be the predicted over- 
lap of the squeezed states at the return time in a perfect 
experiment, then 7 = "f s + j m — ln(a)/i r , where 7 S is 
caused by the laser pulse W, chiefly by photon scatter- 
ing, and 7 m quantifies the decoherence of the cat itself, 
chiefly due to motional effects (electric field noise). We 
see in sets 1-3 behaviour consistent with the expected 
Aa 2 scaling of 7™ [J, 0, 0, EH- We extracted 7 S from 
control experiments at large (200 kHz) detuning, where 
7, n is small. After adjusting for the laser beam power 
used in sets 1-3, the observed value 7 S = 1.7(2) ms _1 
agreed with the photon scattering rate measured in a 
separate test by detecting spin-flips. For data set (3) the 
TDSE gives a — 0.85 so we obtain j m = 3.0(2) ms -1 . 
This is an average over a changing Aa: the coherence 
time T2 = l/2j m — 170(10) /as for this cat state at its 
maximum excursion. 

Numerical simulations showed that the maximum ex- 
cursion reached in phase space, a max , is a function of 
only rj and ao- For ao > 1 a third order polynomial fit 
of the numerical solutions gave 



a max ~ (0.076827a; 3 - 0.45539x 2 



where x = r/ao. Also, the fractional change in the time 
taken for the wavepacket to return to the origin, com- 
pared to LDR behaviour, is linear in the fractional change 
in maximum excursion, i.e. 



1 - R 



tp-t 
to 



0.82 



C*0 "max 

a 



(8) 



We have studied theoretically the forced motion of a 
quantum oscillator subject to a moving periodic poten- 
tial, and experimentally demonstrated large Schrodinger 
cat-like states of spin and motion of a trapped ion, in 
which the return time reveals the non-uniform nature of 
the force, and the inferred motion is such that squeezing 
is confidently expected to be present though not directly 
detected. The coherence time is more than an order of 
magnitude longer than previously reported for this size 
of cat, due to the low heating rate in our trap. The low 
heating rate is a result of the comparatively large size 
of our trap. The coherence time is limited by photon 
scattering and motional heating. The first of these could 
be reduced by detuning the Raman laser further from 
resonance, while the second could be reduced by increas- 
ing the trap size or cooling the electrodes to reduce the 
effect of fluctuating patch potentials on the electrodes. 
The issues studied here are also relevant to quantum in- 
formation experiments where forced motion is used to 
implement 2-qubit quantum logic gates, and high preci- 
sion is essential[l6l[20|,[2ll. 
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(DTO) (W911NF-05-1-0297). 



+1.1352x- 0.011266) /rj (7) 
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